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The Dirac equation with Lorentz violation involves additional coefficients and 
yields a fourth-order polynomial that must be solved to yield the dispersion 
relation. The conventional method of taking the determinant of 4 X 4 matri- 
ces of complex numbers often yields unwieldy dispersion relations. By using 
quaternions, the Dirac equation may be reduced to 2 X 2 form in which the 
structure of the dispersion relations become more transparent. In particular, 
it is found that there are two subsets of Lorentz-violating parameter sets for 
which the dispersion relation is easily solvable. Each subset contains half of the 
parameter space so that all parameters are included. 



1. General introduction 

The concept of searching for small remnant Lorentz- and CPT-violation 
was initially motivated in low energy limits of string theory, but has since 
branched out to include much more general underlying theories of Nature. 1 
Specific coefficients to parameterize the theory were proposed and were 
later generalized slightly to allow for terms violating the gauge invariance 
of the Standard Model. 2 The resulting Dirac equation is modified due to 
the presence of the constant background vector and tensor fields. The first 
general expression for the covariant dispersion relation involving all of the 
coefficients at the same time was initially found using the determinant of a 
4x4 matrix operator. 3 The full expression is rather unwieldy as it leads to 
a fourth-order polynomial in the energy and momentum variables. In this 
work, we use quaternions to reduce the full 4x4 matrix Dirac operator to 
2x2 form. The mathematical properties of the quaternions are then used 
to obtain the dispersion relation in a more tractable form. 4 
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2. Some introduction to quaternions 

The quaternions are generated by four basis elements denoted l,i,j, k, with 
multiplicative properties 

i-j = -j-i = k, (2) 

and cyclic permutations. The quaternions are particularly interesting since 
they satisfy the axioms for a mathematical skew field. A field is an alge- 
braic structure with notions of addition, subtraction, multiplication, and 
division satisfying various properties, the key of which is that there are no 
nontrivial zero divisors. The real numbers provide the simplest example 
with one dimension, they are commutative, and have a natural ordering. 
The complex numbers are a two-dimensional example, but they lose the 
natural ordering property that the real numbers have. The quaternions are 
a four-dimensional field, but they are noncommutative. Octonions are an 
eight-dimensional generalization, but there associativity is also lost making 
them particularly cumbersome to deal with. Frobenius presented a famous 
theorem in 1877 that proves the real, complex, and quaternionic numbers 
are in fact the only finite dimensional, associative division algebras. 

The quaternions may be represented using the Pauli sigma matrices as 

iai — > k, ia 2 — > j, i<r 3 — > i. (3) 

This furnishes an explicit 2x2 representation for the quaternions as 

q = q + qit + q23 +qsk= [ . . (4 

\-q 2 + tq 3 qa - iqij 

The name pure imaginary quaternion is given to a quaternion of the form 
q = q\i + q 2 j + q 3 k where the pure real piece vanishes. Pure imaginary 
quaternions satisfy the following useful relations: 

• ab + ba = —2a • b, 

• ab — ba = — 2 (a x b), 

• abc — cba = 2a ■ (b x c), 

where (a x b) indicates the quaternion that results after applying the con- 
ventional cross product. 

Another useful interpretation of the unit quaternions involves their ac- 
tion as generators of rotations. If q = cos(f5/2) — qsin(8/2), then qaq^ 1 = a' 
where a' = a± cos 9 + (q x a) sin 9 + is the quaternion corresponding to 
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the vector a rotated by an angle 9 about q. The perpendicular and parallel 
components of a are defined with respect to q. Such an interpretation makes 
quaternions useful in areas such as 3D gaming and spacecraft attitude con- 
trols. 



3. Quaternionic form of the Dirac equation 



The gamma matrices in the Dirac representation take the convenient block 
form 



7 = 



a 

-a 



,7° = 



1 

-1 



(5) 



Using quaternions, the conventional free Dirac equation can be expressed 
as 



Po — m 
— po — m 



+ i 



p 

-p 



where 



p = ipjaj = Pai + Pz i + P^ 



(6) 



(7) 



is a pure imaginary quaternion. The lower row gives the spinor solution as 

ip 



p v + m 

and the upper row then reproduces the dispersion relation 

2^2 2 

Pa - V = m . 



(8) 



(9) 



Note that <j> and x arc two-component complex spinors that serve as a 
module on which the quaternions act. This is in contrast to other approaches 
to using quaternionic valued wave functions as has been attempted several 
times in the literature. 



4. Quaternionic form of the perturbed Dirac equation 

The Dirac equation with Lorentz violation takes the form 

(r% - M)ip = o, (10) 

where 

r" = 7 " + c^7m + <^757m + e" + if v J 5 + (11) 
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and 



M = m + + b^-f + \B^ V . 



(12) 



The coefficients a^, 6 M , etc... are Lorentz- violating constant background 
fields. 

Some field redefinitions can be used to simplify the model. 6 The first 7 
is to fix r° — » 7° using the transformation ip = Ax, A = (7°r ) -1 / 2 . 
This makes the hamiltonian hermitian and is often a necessary first step 
in properly interpreting the Lorentz- violating physical effects of any given 
experiment. The second involves a linear transformation on the momentum 
and mass space to remove o M , , and e M terms 



P =p-a-c p , 

p' = (1 + c 00 )po -a -C!-p, 

m! = m — e ■ p. 



(13) 
(14) 
(15) 



Since the transformation is linear, it is simple to invert at the end of the 
calculation. To avoid cumbersome notation, the primes are dropped in the 
rest of the calculations. Finally, the / M term can be removed using a more 
complicated redefinition. 8 

The perturbed Dirac equation can then be put into quaternionic form 



p a - m a - 
-a + & ~Po 



a 

- m 



+ i 



6 ~? 
p 6 



(16) 



The parameters in the Dirac operator with hats are pure imaginary quater- 
nions, while the unhatted objects are real. Their relation to the original 
Lorentz- violating coefficients are given by 



a = b° + di ■ p, 
e = b + d p + (g - h), 

with the notation 

Mi 



a = H — G, 

6 = -b-d p + (g - h), 



(17) 
(18) 



d\ = d v 
K l = l/2e tlk W k , 



C - 9° ij Pj 
9 



l/2e ijk g ikl Pl . 



Substituting the first row equation of Eq. (16) into the second row yields 



(a - a + ip)(x - ie)(ao + a + ip) + r(—y + iS) 



(19) 
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where x = p° — m, y = p° + m, and r = x 2 — |e| 2 . This equation reduces to 
the form 

(qi + m)i = 0, (20) 

where qi is real and q 2 is a pure-imaginary quaternion. This equation yields 
the dispersion relation from the eigenvalue condition \q 2 \ = q\. A fourth- 
order polynomial in the energy p° results in the form 

A) + PiPo + fa(po) 2 + (Po) 4 = 0. (21) 

The linear term fi\ vanishes for special parameter choices leading to a simple 
factorization of the dispersion relation. 

Case 1 involves setting e = S, a = 0. This is equivalent to setting 
b = dp = H = G = 0, while leaving the other parameters arbitrary. The 
resulting dispersion relation takes the form 

p 2 = f + m 2 + a 2 + S 2 ± 2 A /£>i(p), (22) 

where 

£>i(p) - (5 'x pf + (a P ~ m5f . (23) 

Case 2 involves setting e = —5, ao = 0. This is equivalent to setting 
h = g = bo = d\ = 0, while leaving the other parameters arbitrary. The 
resulting dispersion relation takes the form 

p\ = p 2 + m 2 + a 2 + S 2 ± 2y/D 2 (p), (24) 

where 

D 2 {p) = (axp-mS) 2 + (5-p) 2 + {a-5f. (25) 
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